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BACKCASTING FROM GRADUATION TARGETS TO 
REQUIRED ENROLLMENTS USING THE GENERALIZED ■ • 

INVERSE OF THE TRANSFORMATION MATRIX 

The general solution to the problem of calculating school 
-.enrollments at time (t-1) necessary to meet predetermined graduation or 
manpower targets at time t is derived anc^ examined. It/s based on 
the generalized inverse of the transformatioQ^ia^i^^hat describes 
the pupil flow through the school. The-author has applied this technique 
[Jaoua et al, 1974] but at the time therV^/as no solid theoretical 
basis for it. Tills paper provides one. 

The field of general ized ^inverses xbntains a confusion of definitions 
and names including conditional inverse, pseudoinverse, and weak 
generalized inverse. The one used here is the stong or unique generalized 
inverse defined by the four equations below (1 to 4). For many but not 
all of -the findings in this paper a weaker nonunique inverse meeting 
fewer than 4 of the conditions would have sufficed. To achieve sim- 
plicity of presentation and\because of the availability of computer 
software we use only the unique one. 

Symbols and conventions are now outlined. "A" represents a 
general matrix which transforms enrollments into graduations. "fc M 
1s Its generalized Inverse. "B" and "C" are two matrtces whose product 
1s A (A =BC) and whose Inner dimensions are equal to the rank of all 
three, "x" 1s a vector of enrollments and "V a vector of graduates , 
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each disaggregated by discipline, program or labor farce category. 
V arid V are scalars representing vector or matrix dimensions, 
"r" 'always represents the rank of a matrix and sometimes a dimension. 
Lambdas ("A") are eigenvalues; "T is the identity matrix; - and "h" 
is a vector of arbitrary variables .jefith elements h,, h 9 > etc. "Iff 11 
means if and only if.^ The word "solution" is, used rather loosely in 
that it sometimes refers to an exact genuine solution and sometimes to 
an approximation. The precise meaning is clfear from the context. 



a^k( 



The Backcasting Problem . - ■ • 

Numerous authors have modeled student flow through school systems 
[See Johnstone, 1974; and Schiefelbein and Davis, ,1974 for reviews] f 
using rtiatrix transformations of the type Ax = y v/here is the prd- 
portion^of persons in grade>or institution j (often called "state") 
fc during the initial time period who move onto state 1 for a later time 
period. The vectors x and .y represent the number of persons in the 
various states during the initial and later periods respectively. If 
the transformation matrix A is square, and equal i's and represent 
the same states then the transformation can be compounded over any number 
of periods by raising A to higher powers. In the case of square A this- 
is similar to a Markov process except that it is frequently seen as 
deterministic or prescribed rather than stochastic. Furthermore, drop- 
outs are often treated as a residual in which case [l^A t [T]* and 
populations are not conserved. / 

Frequently educational planners are given this problem in 
reverse. They are provided with a set of manpower /or graduation targets 
to be met in some future period^ They must then calculate the enrollments 
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necessary to bring about the required number of graduations. If the 
transformation matrix A has equivalent row and column states and is 
square and nonsingular then the solution is obtained by premul t1 plying \ 
the graduation target vector y by the Inverse df A (x = A^y)* It 1s 
argued below that A 1s usually notlsquare and that 1f 1t 1s made square 
by Including all states 1n both the row and column spaces 1t will, 1n 

i 

general , be singular. 

Any square transformation matrix that accurately represents 
progress through a school or school system is nllpotent and therefore 
singular. Everyone who enters a school eventual l^teaves either through 
graduation, dropout, or death. Thus some power of the matrix A which 
transforms entrants Into graduates must map every entrant vector Into 
a zero graduation vector. The only such matrix 1s the -zero matrix. , 

If A n = 0*then A i$ nllpotent which 1s sufficient for singularity. 

/ * 

The addition of an absorbing state to the matrix to collect all 

the graduates prevents nil potency but does not guarantee 1nvert1b1l1ty. 
In general, If the. system contains division points where students flow 
from'one stream equally Into two or more streams the transformation 
matrix 1s singular. 

i The problem Qf backcastlng from a vector of graduatlpn targets 
to the prerequisite^ enrollments requires a. more general formulation ' 
than that provided by the standard Inverse for a nonsingular square 
matrix. This paper shows that the problem can be solved using the 
cqncept of a generalized inverse matrix [Ben-Israel and Grevllle, 1974; 
Bouillon and Odell, 1971; and Rao and M1tra, 1971] with properties 
expressed by £qua,t1ons (1) to (4). A 1s any matrix and 6 1s Its gen- 
eralized Inverse. * 
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AGA = A (0 

GAG = G (2) 

AG = (AG)* i (3) 

GA = (GA)* * (4) 

SSI • *~ 

The generalized inverse G exists : for all A and 1s unique /.Grayblll., 

* -l ' 

1969, p. 97.]. Substitution of A for G in the above shows that the 

definition of G simplifies to the standard Inverse for 1nvert1ble A. 
For rectangular m x n matrices, G has dimensions n x m. The Inverse of 
a' zero matrix 1s the zero matrix- with reversed dimensions. 

The following two sections develop a- "generalized Inverse" with 
the properties needed to solve the backcastlng problem. This 1s then- 
shown to be the matrix that satisfies conditions (1) to (4) above. 

7 
> 
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The Least Squares Criterion * . 

7i 1 

' ' The starting point 1$ the matrix equation Ax ■ y with A, a 

given m x n transformation fnatrlx of rank r; and y, a vector*of*7jradua- 

tlon targets. The solut1oij( set which contains values of the enrollment 

vector x may be empty. If^ 1t 1s, an approximation to x 1s desired. 

The approximation used 1s least squares, the sqlutlon that( 
minimizes the squared discrepancy between the actual y and the y pre- 
dicted by x. We wish to "Select an x that\m1n1m1zes 

(y-Ax) t (y-Ax) ' ( ' (5) 



|j (y-Ax) t (y-Ax) = ^A^y-Ax) = 0 



(6^ 



A*Ax -AV , (7) 
A deus ex machlna 1s required at this point because A may not be ,1nvert1ble 
This 1s provided by thejsingular value decomposition of A [Noble, 1969, 



p. 338]. A 1s decomposed Into two .matrices B and C of order m x r and 

r x n respectively such that A = BC. All tftree matrices are of rank 

"■'.•.„• > 
t- t ' 

r. For this reason (B B) and (CC ^are 1nvert1ble. Substitution of 

BC for A gives: 

(BC) t (BC)x - (BC)*y (8) 
Premultiplylng by C and regrouping gives: 

m 

(CC t )(B t B)Cx = (£C t )B t y (9) - 

Because of the Invertlfclllty of the leftmost terms the expression reduces 
to: 

Cx = (B^rVy (10) 
Equation (10) gives the least squares solution for x as far 
as v/e can take it. It will always exist but may not be unique. There 
may^be^an Infinite set of x's that reduces the squared discrepancy to 
2ero or to some single positive value, ✓ 

The Minimum Norm Criterion 

To obtain a unique solution another criterion^ 1s needed. We 
choose, to minimize the Euqlldean norm of x, using the techniques of 
Lagrange multipliers. 0 ■ ■ 

^ , L - x*x + 2X*(y* - Cx) (11) 

where y* « (B t B)" 1 B t y \ (12) 

The symbol y* defined by equation (12) 1s used for convenience. The 
constraint 1n parentheses 1n (11) 1s just equation (10). The , term 2X 
1s a row vector of Lagrange multipliers with the scalar , 2, also present 
for convenience. Partial differentiation yields equations 413) and (14). 



ft- = 2x - 2CH = 0 iv (13) 

ft. 2(y* - Cx) = /) (14) 

From Equations (13) and 04) we obtain (15) and (16) respectively. 

x ■ C*X - , . (15) 

y* -iCx ,J (16) 

Substituting (15) Into (lw gives: 

y* = CC** | . (17) 

X -'(CcVV : . • < 18 ) 

Substituting (18) Into (15) gives . \ 

. x = cVrV'' . . , *(19) 

1 Finally, we substitute (12> into (19) to get rid of y* and give the least 
'"** i • 

squares minlmunf norm solution. 

x = C t (GC t )" J (B t ^)" 1 B t y ' (20) 

♦ Equation (20) Is! now used to define G. 

. Q C^CC*)*- 1 ^)"^* V • (21) 
' Substitution of BC for A and the fclgjvt^lde of (21) into equations 
(1) to (4) will convince the reader that' the above expression satisfies 
the definition of the unique general 1zid^ Inverse. 

S ome Generalized Inverse Properties 

■ ■ WW •-MT'I" ' " * ' ' — — — — 

A number of Important results which afife used later follow from 

TV / 

the definition of 6 (21) and the fact that/A - BC. These are listed 

,t j 

below:., " i 

>V AG - B(B t B)" 1 B t i ' (22) 

* GA - C^CcVV I • ' (23) 



■'I ,„„x2 



V 



(GA) - GA > 1 , (24) 
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Note that AG and GA are m x m and n x n matrices, respectively. Equation (24), 
which. can be derived by squaring (23), shows that GA is an idempotent. ■'. > 4 
matrix. Premult1pl1cat1on of (24) by (GA)" 1 demonstrates that the only 
nonsingular idempotent matrix 1s the Identity matrix. Squaring of ^ 
(I-6A) followed by the application of (24) can be used to prove that 
(I-GA) 1s also idempotent* an important result for the next section. 
Similarly, (AG) and (I-AG) are idempotent^ The expressions for GA and . J 
AG can be used to show that GA f AG f 4 (I-GA), and (I-AG) are all symmetric. 

The expression for G (21) gives us the least squares minimum 
norm solution to the equation Ax = y as follows: 

^ x.- Gy' ' <25) 

# But what is really needed 1s a solution to equation (10), that tt 

is, the set of all x for which the squared discrepancy, ((y-Ax)*(y-Ax)) f 

t reaches its minimum value. If the minimum squared discrepancy 1s zero, 
the solutions are^exact; otherwise they are approximate. Either may 
be unique or Indeterminate. The general least squares solution is 
presented 1n the following section followed by a proof that It exactly 
represents the complete solution 'to equation (10). 




The General Least Squares Solution 

The complete least squares solution to x is as follows: 
x » Gy + (I - GA)h / . (26) 

where h is an arbitrary n x 1 column vector. This is proven in two* 
steps. Firstly, it 1s proven that (26) satisfies (10) and secondly, 
that for every x satisfying (10) there is at least one h such that x 
also satisfies (26). 

1 u 



The first step begins with -the substitution of (26) into (10). 

"c(6y -6A)h) « (B*B)" 1 B t y.. - ° ( 27 ) 

Substituting (23) for GA and (21) for 6 gives the following expression 

for the" left hand side 

CjC t (QC t )' 1 (B t B)' 1 B. t y + (I - C t (CC t )" 1 C) h) (28) 

t t -1 

Removing the outer parentheses and removing, CC (CC ) 1n two places 
results 1n considerable simplification. 

(B^rVy + (CI - (?) h # (29) 

The right .term cl.early disappears ^leaving (29) Identical to the right 
hand side of (27). Thus our expression for x, (26), 1s a solution to 
(10). .We now prove that (26) 1s the complete solution.' 

Starting with equation (10) we show that every x that 1s a 
solution 1s also a solution to (26) for some h. Premult1ply1ng (10) 
by GB gives: 

.* i * " '• * 
GBCx 1 = GB(B l B)"Vy * (30) 

Substituting A for BC In the left and equation (22) 1n the right gives: 

GAx = (GAG)y (31) 

Substituting (2) Into the right and rearranging yields: 

0 = Gy - GAx (32) 

x = Gy + x - GAx (33) 

"x a Gy + (I - GA)x (34) 

Thus any x substituted for h satisfies equations (26) and (10). This 

does not Imply that x can be generated only by Itself. The h corresponding 

to any given x 1s usually hot unique. ' 

In the formula x ■ Gy + (I - GA ) h , the first term on the right 

1s the solution corresponding to h, s 0. It Is also the minimum norm 

solution as shown in (25). The second tern; provides for all other 
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solutions. Interpreting this as a solution to a system of linear equations, 
we see that the first term 1s a particular solution (exact or approximate) 
of the Inhomogeneous system and the second termMs the general solution 
of the homogeneous system. Attention 1s now concentrated on the second 
term 1n order to determine whether the solution 1s unique, and 1f not, 
how many elements there are 1n the arbitrary vector, h. 

The Indeterminacy of the Solution 

. x In this section 1t is first shown that the rjnk of GA 1s equal 
'to the sum of Its dl^&nal elements, called Its trace. The proof depends V- 
.only on GA's symmetry and Idempotency so that we can also apply the 
finding to I and (I-GA). Thfs is followed by three simple proofs, that 
rank(GA) ° rank(A), that rank(I-GA) ° n - rank(GA) and finally that 
rank(I-GA) « nulHty(A). We^ conclude that the columns of^(I-GA) form 
a basis for the null space of A and that the rank of t (I-GA) equals the 
minimum number of elements 1n the arbitrary vector h. 

Since GA 1s real, square (n x n), and symmetric 1t 1s 4 Herm1t1an, 
from which a number of Important results follow. The eigenvectors of 

*# 

GA are all linearly independent and, 1n fact, are orthogonal [Noble, p. 321]. 

4 » 

f 

GA can be reduced to diagonal form by a transformation that is both a 
similarity and a unitary transformation. 

D = Q" 1 (GA)Q - ' (-35) 

or D ■ Q t (GA)Q . ' (36) 

such that, 

Q t f Q' 1 ^ ' f . (37) 

Q 1s an orthonormal matrix and D 1s an n x n. matrix containing the\^ 
eigenvalues of GA 1n the diagonal and zeros elsewhere ENoble, 1969, p. 3^8]. 

■ "V 

1 ^ 
1 *~ f 



)'. 



By prertiulti plying the equation for eigenvalues, (GA - 'x I) v. *•-()•, 
ty GA arrcJ substitutinq GA for (GA) the deader # can derive the equation 
(1 - x) (GA)y - 0 and conclude that either X - 1 qr'GAv = 0, in which 
case X = 0./ Thus the. eigenvalues are all zeros and ones. 

Consider the relationship between the' ranks and eigenvalues of 

QA ^nd D. The rank of D is equal to the rank of ;GA/because multiplication 

' -1 
of a maiV'ix (6A) by nonsingular matrices (Q and Q) does not change its 

rank [Nlble, 1969, p. 138], The rows of the diagbrial matrix D can be 

interchanged to put the identity matrix into the upper left corner and 

zefbs lierywhere else. From this it becomes clear that the rank of D 

is'equ$ to the sum of its eigenvalues. The same conclusion holds for 

GA since D and GAhave the same eigenvalues. 

. •' ; /■ ■' V . v -' >, -,.f > ■ ' ' 

This brings us to an important finding which applies to all ■ ■ 

. " . " : , / . ■ ' . v ' : • ■ ; : • 

symmetric idempofent matrices Because the sum of the eigenval ues of 

ajyu^uare matrix is equal to its trace [Graybill, 1-969, p. 223], the 

rank or GA is equal toWts trace,. This provides a simple computational 

technique for determining the rank of GA and (I-GA), both of whiclf^are 

symmetric and idempotent. . </ / \ 

* ■ * . ■ • .. ■'<•»<► * ■ . " 

Since the tr^ce of a matrix difference is equal to the difference 
of the traces, the rank of (IrGA) is equal to v rank(I) > raijk(6A) or 
n - Tank(GA). 

■7 It is well known that for any G and A confprmable for multiplicatio 
rank(GA) <_ rank(A) [Noble, 1969, p: 139]. Applying this same principle 
to equation (1} grouped as A(GA) s , A,< we caji, conclude that rank(GA) >_ rank( 
Therefore, rank(GA) = rdnk(A)..> By similar reasoning it can be shown that 

1 . t V \ 

all of A, G, GA, and AG'hatfe'the same rarfic, as do B and C. x 

We can now combine these findings to reach the final conclusion 



Of this section. The general solution contains a matrix (I-GA) whose 
rank is equal to its trace. Furthermore, its rank is equal to n- rank(GA)- 
or n - rank(A). But n - rank(A) = nullity(A)' by Sylvester's Law of Nullity > 
[G1llett, 1970, p. 154]. • Therefore, the trace and rank of (I-GA) are \ 
equal to the nullity of A. By premultiplylng (I-GA) h" by A and substitute 
A for AGA (2) we see that all vectors of the form (I-GA)h«are mapped 
Into the zero vector by A. Thus, tbejftre in the null space of A. 
Because (I-GA) has the rank of the null space of A, (n - rank(A))*, 1t . 
must have exactly (n - rank(A)) linearly independent columns. The 
col-umns of (I-GA) therefore form a ba$is for the null spaqe of A. 

Unfortunately, this is not a minimal basis for the null space 
of A; 1t contains n vectors 1 but only (n- - rank(A)) linearly Independent 
ones. This form of (I-GA) does not ^provide much insight .into the structure 

' 1 ■ ' . . . ' 

of A's null space, a practical and computational problem which is discussed 
In the section on computation. 

The Approximateness of the Sol ution & - . 

■ > - . . % r . * * ■ 

This section examines the circumstances underiwhich the equation. 
Ax=y yields exact and* approximate values for x. Firstly,-, an expression 
for the squared discrepancy that depends only on y and A is derived. 

This is used to define three possible outcomes for a given A and the set 

• \ ., • ■ ' 

Of all y. i 

. .**••• i ■■§>■ 

The squared discrepancy 1s ((y-Ax) (y-Ax)) which 1s minimal* 

when x = Gy + (I-GA)h. Substituting the' expression for x 1n the expression 

for the squared discrepancy gives: " ; 

^ !> (y-A(/y + (I-GA)h^(y-A(Gy + (I-GA)h)) * ' - , (38) 

This expression can be simplified by substituting zero for A(I*-GA). 



• , (y < rAgy) t .(y-AGy) • ' ; • (39) 

yV - y t (A6y) < - (AGy)^ + (AGy) t AGy (40) 
This can be° simplified by looting that each <term is «a scalar and, therefore, 
equal to its, transpose and that AG is symmetric and idempotent.' 

V(I-AG)y ' x • ' (41) 

Equation (41) 1s the expression we seek, the minimum squared discrepancy 
for a given A and y. Exact solution(s) exist wh~en (41) is equal to zero. 

To examine the conditions under which (41) is equal to zero 1t 
is necessary to use the fact that a nonnegative quadratic form is equal 
to zero if and only if -the vector valued expression produced by removing 
the prepositioned Vector, from the quadratic form' is zero. 

v The proof of) this involves the diagonal ization of (I-AG)^and the, 

/ • t t 
transformation of y and y Into z and z by the orthonormal matrix of 

eigenvectors. ■ j- ' 

y t (I-AG.)y = yVW - (te) 

■ t We define z as Qy. Since Q .is orthonormal z*'= y*Q _1 . The express ion J* 

becomes s 

z*Dz -. ' (43) 
where D is the diagonal matrix of eigenvalues of X I -AG ) . This expression 
is a scalar and can be written as follows: , 

z^z = ZjUjZj) + z 2 (X 2 z 2 ) + ... ' (44). 
where the subscripted z's represent elements of the vector z. Clearly 
(Zj(X]Zj)) = 0 if and, only if ■ 0. Since the eigenvalues are all 
nonnegative the expression (ZjXjZ.) must also be nonnegative. Thus, 
the quadratic form can equal zero only if every term on . the right hand 
side of (44) is. zero, which means that every x,z. = 0. This is the key 



to the proof. It means that: 

'"ySj^DQy = 0 Iff DQy = 0 \ 
-1 

Since Q is nortsingular, 

PQy = 0 iff Q -1 DQy:= 0 
Combining (45) and (46) and substituting .(NAG) for Q 'DQ completes the t . 
proof. ,* * 

y t (I-AG)y = 0' iff (I-AG)y = 0 ■ 

We can now proceed to examine the conditions under which exact, 
solutions exist.- There is an exact solution formal 1 y such that 
S y .= AGh * ' (48) 

The" first stage of the/proof verifies that any y satisfying (48) results . * 
In a zero value for (41), the squared discrepancy. Substituting (48) Into 
(41 ) gives: \ . r 

(AGh) t (NAG)AGh (49) 
Multiplying (49) through by the post positioned AG and application of'(T) 
results 1n a zero, value for the expression. The second stage of the proof , 
shows that art h exists for every y that results in a zero value of (41). ^ 
This time the expression (I-AG)y from (47) 1s- set to, zero, rather than the 
corresponding quadratic form. Manipulation of this expression quickly / 
shows that: 

y - AGy ' ' •* ' (50 

This completel the proof. Since AG 1s generally singular, y 1s not the 

% . ■ 

only value of h that will generate y. 

. The results for solution approxlmateness can t£ summarized 1nto^ 
three cases, the first of which 1s trivial. 

If AG*0 the only y yielding an* exadt solution 1s the zero 
vector. Inspection of (l) shows that AG=0 1f and only 1f A=0. 
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From equation' (41) 1t is clear that if AG = I, there fs an 

exact'solution for all y. Since AG is m x ra /nd idempotent this can ■ 

. occur only jf the rank of A 1s m. 

The intermediate cases are those for which the rank of AG 1s 
V / / 

greater than zero but less than m. These y;eld exact. solutions for 

/ * • • . 

those y that satisfy (50) and approximate .solutions for all other values 



of y. 



Computation 

This section contains/an example of an enrollment backcast and 
some hints on computation. 

The example was s/lved with a program of 50 Fortran statements / 
that Input the transformation/matrix (A), the target vector (y), and then 
called widely availab/e subroutines to do the rest. The subroutines, 

Including one that /alcul/tes the unique generalized inverse of a matrix 

./.<./' 

by means of singular value decomposition. [Noble, 1969, p. 335-40], are a 
part of the International Mathematical and Statistical Library [IMSL, 197.4]. 
-This library a^so contains subroutines to do the other matrix and vector 
operations necessary to calculate and output the complete solution. - 



*** Place Figures 1 anck2 about here *** 

The problem and solution are contained 1n .Figures 1 and 2, . 
respectively. The problem consists of one transformation matrix (A) o 
two graduation target vectors, y,, and y 2 , which generate two least 



an 
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solutions, x 1 and x 2 , respectively. The general solution to the 



'homogeneous problem 1s the same for all y, given A. 
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The information in Figure 2 permits us to analyze the solutions. 

The matrix (I-GA) is 4 x 4 with a trace of 2. Its rank is, therefore, 

2 and from this GA has a rank of 4-2*2. V A and AG are of rank 2. (I -AG) 

is 5 x 5 and has a rank of 5-2=3. We conclude that there will be an 

infinite number of least squares solutions, some exact and some approximate 

t * 

depending on y. The expressions for y (I-AG)y show that y-j has exact 

solutions and y 2 has approximate solutions. 

The matrix (I-tSA), which is a n&djjndant basis of the null space 

of A, can be converted to a more convenient form by removing the unnecessary 
-vectors-. A systematic way to^do this is to reduce 1t to row echelon form 
Tfe lso k nown as Herm^e normal- form) ^d transpose it. -After discardin g 

the zero columns this gives a basis of (n-r) column vectors with the "* ': 

Identity matrix 1n the top (n r r») rows as shown at the bot'tom of Figure 2. 
Two undesirable types of enrollment vectors are often generated. 

The least serious and most common contains nonipteger ervrol 1 men ts that «i 

can be rounded. Negative enrollments" mjty occur inVirher the leas£norm - 

* 

or £he general solution. „ 

. ■ ' s " ' r 9 ' ? 

Th£. numbSr of enrollment vectors 1n any solution wf 1.1 be finite 

if only vectors with nonne^atlye integer enrollments with a finite maximum 

sum are considered. * /• 

* ^ * \ « 

5* '\ 

* \ ft 

* ■ t » . 

Conclusion ^ " ' 



The enrollment backcastlng problem Ax=y has been so Wed t The 
complete exact or approximate solution 1s x = Gy + (I-GA)h where G is the 
generalized Inverse of A and h 1s an arbitrary vector whose mlnlmiiiTv 
dimension 1s equaT to the nullity of A. The solution 1s- exact 1f 
y t (I-AG)y =0; otherwise, 1t 1s a least squares approximation. Figure 3 
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FIGURE 1 
THE EXEMPLARY PROBLEM 



A = 



.04* 


• SCO 


.2b4 


• 13S 


•210 


• 050 




.162 


.44U 




.195 




.1 10 


•own 


•313 


.197 


;20B 


• 100 


• 138 


.17& 



?32 . 0 
1 14.0 

204.0 
140.0 



1%.0 
124.0 
10.0 
^5S.0 
334.0 
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FIGURE 2 

THE SOLUTION TO THE EXEMPLARY PROBLEM 



U~l. ?*!■ 1034491 
1H.5!, ^24 13733 



x 2 - 



?u?40 



7. P 150 

4: £.2*737583. . 
*n*. M !SC^753 
I ?6. 358 7207460 



y 1 t (NAG)y 1 - §; 



y 2 (I-A6)y 2 * [5s344.rj 



(I-GA)h 



.1189555 172 
. IS 183*1828 
.077:662069 
.4524137931 



.1810344828 
.31S9SS51F2 
.4224137931 

.0775862069 



-.0775862069 
-.4224137931 
.6810344928 
IS 10344323 



-.422413/$31 
••07.75£i§t069 
•. 1010^S#2g 




(I-GA)h 
reduced 



1 . 0090000000 
0. 

. 75OS0OOH00 
• 1 . 7300000000 



0. 

1 . 0000000000 

■ 1 . 7^00000000 

. 7500000000 
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ANALYSIS OF THE GENERAL SOLUTION 
BY MEANS' OP RANK AND ORDER 
INFORMATION 



Start here. 







A 


is m x n 


and, 


of raftik r. 




n>m=r , AG=I 
There are an infinite 
number of exact solu- 
tions fot every y 



T<m and r<n 
There are an infinite number of 
solutions which are exact 1J£ y=AGy 
and approximate if y/AGy. 




/m=n^, G=A 
Thete is- unique 

exact x for every y. 

— 



n£>n=r, GA«I 
Hiere is a unique solution 
which is exact if y^AGy and 
approximate if y/AGy. 



NOTE: "G" is the unique" generalized^ inverse of A. 



2.1 



